Is the ground state of Yang-Mills theory Coulombic? 
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We study trial states modelling the heavy quark-antiquark ground state in SU(2) Yang-Mills 
theory. A state describing the flux tube between quarks as a thin string of glue is found to be a 
poor description of the continuum ground state; the infinitesimal thickness of the string leads to UV 
artifacts which suppress the overlap with the ground state. Contrastingly, a state which surrounds 
the quarks with non-abelian Coulomb fields is found to have a good overlap with the ground state 
for all charge separations. In fact, the overlap increases as the lattice regulator is removed. This 
opens up the possibility that the Coulomb state is the true ground state in the continuum limit. 



I. INTRODUCTION 

Explaining colour confinement remains the outstand- 
ing problem in the theory of strong interactions. The con- 
finement mechanism is typically pictured, in the mesonic 
sector for example, as resulting from a change in the dis- 
tribution of the glue around two quarks as they are moved 
apart. It is expected that the glue forms a flux tube, 
the exact geometry of which is unknown, connecting the 
charges and resisting their separation. This idea has been 
supported by numerous lattice simulations which have 
revealed that at large separations the ground state po- 
tential increases linearly with the distance r between the 
quarks (see for example 

In this paper we will model the confining flux tube 
by generating gauge invariant states | \& ) which de- 
scribe possible configurations of the gluonic fields around 
charges, and testing how well they overlap with the un- 
known ground state in the mesonic sector, |fi). We will 
build and test various proposals for such a state, using 
a combination of analytic and lattice techniques. Our 
focus in this paper will be the cutoff dependence of the 
models, as imposed by the short distance lattice regula- 
tor, and the behaviour of our models in the continuum 
limit. 

Central to the construction of our trial states is the 
need for gauge invariance: a naive construction of a 
mesonic state, where we act on the empty vacuum | ) 
with fermionic creation operators, fails to produce a 
physical state - it is not gauge invariant. The origin 
of the problem lies with the Lagrangian fermions, which 
themselves are not gauge invariant and do not immedi- 
ately describe physical observables. In fact, no such local 
operator can create physical charges in a gauge theory Q . 

The basis of this observation is that charges always ap- 
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pear with an associated (chromo-) electromagnetic field. 
The result is that the physical electron, for example, is 
a composite object comprised of a U(l) fermion together 
with its surrounding electromagnetic fields. This mani- 
festly non-local object is gauge invariant and describes an 
observable charge [3, 0, U Q ■ Quantum mechanically, we 
see that to construct a charged state, or a multi-fermion 
state with zero overall charge, the fermion creation oper- 
ators must be accompanied by a configuration of gauge 
boson creation operators, which generate the necessary 
chromo - electromagnetic fields. We refer to constructing 
these fields as 'dressing' the charges. Although gauge in- 
variance imposes restrictions on such a state, there is still 
a huge choice in the form of the dressing. 

It transpires that any gauge fixing condition may be 
used to define a dressing and therefore a gauge invariant 
state @. Additional conditions on the state can further 
refine the class of gauge fixing used 0]. In particular, the 
Coulomb gauge is particularly relevant to the description 
of static charges and this can be clearly understood per- 
turbatively [5[ . The extent to which this is also the case 
non-perturbatively is the focus of this work. 

In this paper we concentrate on two model states re- 
lated to the Coulomb and axial gauge fixing conditions. 
The 'Coulomb state', |<E>), describes two individually 
gauge invariant charges each surrounded by non-abelian 
Coulombic fields. The 'axial state', which we denote | % ) 
throughout, describes an overall gauge invariant meson 
formed by linking two fermions by a string of glue. We 
will formalise these descriptions below. Here we note 
two properties of these states: first, although our nu- 
merical calculations are in SU(2), our analytic construc- 
tions apply to SU(A C ); secondly, if we represent the states 
as wavefunctionals of (the matter fields and) the gauge 
field A, then the dressing becomes unity when A satis- 
fies the gauge condition defining the dressing. For ex- 
ample, the Coulomb dressing may be written as a func- 
tional of A which is always equal to 1 when d{Ai — 0. 
Therefore, if we work in the gauge used to define the 
dressing, then dressed fermions are just the usual bare 
fermions. This observation provides a link between lat- 
tice studies of manifestly gauge invariant states and ana- 
lytical techniques such as the Hamiltonian functional ap- 
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proach based upon Coulomb gauge @, 0, H, H, E3, [0, EH , 
and investi gat ions of the non-abelian Coulomb potential 

Hi Q EH 0H. 

Once we have constructed a trial state "J ) we compare 
it with the unknown ground state through their (mod- 
squared) overlap | ^f)\ 2 . This overlap may be calcu- 
lated on the lattice through the persistence amplitude 
( ^ \ e~ HT \ ^ }. Since the large Euclidean time limit is a 
ground state projector, we have 

(*|e^ T |*) ^e- v ^ T \(n\^)\ 2 , T»l, (1) 

where the ground state energy in the q-q sector gives 
the well-known confining potential, V(r). This holds 
provided that the overlap between the trial and ground 
states is non-zero, which has recently been verified for 
our trial states [1 9j : the large time persistence ampli- 
tudes for both the Coulomb and axial states yield the 
confining potential V(r), confirming a non-zero overlap 
with the ground state. 

In [2(| we studied the overlap for our trial states in 
U(l) and SU(2) Higgs theories, in both the confining 
and deconfined phases. In the U(l) theory with heavy 
charges exact calculations can be performed in the de- 
confined phase. There we found that the ground state 
is exactly the Coulomb state $ ) , describing two indi- 
vidually gauge invariant charges with the familiar elec- 
tromagnetic Coulombic fields. We found that the axial 
state | x ) is an infinitely excited state, giving, in the time 
slice where it is prepared, a confining potential between 
the charges. The state is unstable and decays in time 
to the ground state. Contrastingly, we found that in the 
confining phase of U(I) the axial state had the larger 
overlap with the ground state. The physics in this phase 
therefore appears to be better described by the confining 
potential generated by a very thin string. 

In the Higgs phase of SU(2)-Higgs theory, we found 
that the Coulomb state had the better overlap with the 
true ground state. However, contrary to what may have 
been expected, the Coulomb state also provided a better 
description of the ground state in the confining phase 
than the axial state. This implies that the SU(2) flux 
tube is significantly thicker than the confining string in 
U(l) theory. 

Our previous simulations were performed at a fixed lat- 
tice spacing. In this paper we will investigate the depen- 
dence of our results on the cutoff and so study the contin- 
uum limit. We begin in Section [TT] with the axial state, 
which describes the flux tube as a string of glue, and 
comment on some recent results in the literature. In Sec- 
tion lllll we compare our results with those for a thickened 
string constructed from smeared links on the lattice. We 
will see that the overlap between each of these states with 
the ground state decreases as we approach the continuum 
limit, indicating that they have little to do with ground 
state physics in the continuum. In Section HVl we turn to 
the Coulomb state, which describes a charge-anticharge 
pair surrounded by non-abelian Coulomb fields. We find 
very different behaviour as the lattice regulator is re- 



moved: the overlap of the Coulomb and ground states 
now increases as we approach the continuum. Finally, in 
Section [V] we present our conclusions. 



II. THE AXIAL STATE 



A. Construction and properties 



The axial state is 



X) = g(a; 2 )Pexp 



3(*i)|0>, (2) 



with q(x) the heavy fermion source. This state describes 
a gauge invariant, charge neutral meson formed by link- 
ing the matter fields by a gluonic string. We take the 
string to lie on the straight line connecting x\ and X2 1 
and so r = |a;2 — £Ci| is the separation of our fermions. 
For the axial state, calculating the persistence amplitude 
over time T corresponds to calculating an ordinary rect- 
angular (and unsmeared) Wilson loop of spatial extent 
r and temporal extent T |19| . Before presenting the nu- 
merical results it will be useful to outline the behaviour 
of the analogous state in QED with heavy charges. There 
the ground state is known exactly, and its overlap with 
the axial state is zero [20| . The origin of this vanish- 
ing overlap can be traced back to short distance effects, 
as the flux is trapped on an infinitely thin string, even 
though the axial wavefunctional contains no UV (nor IR) 
divergence. Including a momentum cutoff A, the overlap 
is found to be, up to terms finite as A — > oo in the expo- 
nent, 

|(fi| x )| 2 = (rA) 4< */ ff exp[-arA + ...] . (3) 

Here a = e 2 /(47r). The overlap is a function of rA, dom- 
inated by exponential decay for large values of the cutoff 
which probe the infinitesimal extent of the string in the 
two directions transverse to X2 — x\. As the cutoff is 
removed this causes the overlap with the ground state to 
vanish. 

We now describe our calculation of the unsmeared 
SU(2) Wilson loop. We used a standard heat bath algo- 
rithm combined with microcanonical reflections to gener- 
ate an ensemble of 1000 configurations which were used 
to measure the persistence amplitudes. Defining a super- 
sweep to be a combination of 3 heat bath update sweeps 
followed by 7 microcanonical reflections to enhance the 
ergodicity of the algorithm, we allowed 250 supersweeps 
to reach thermal equilibrium. We found no dependence 
of our results on the initial configuration chosen for the 
lattice. All measurements were taken on a 20 4 lattice, of 
varying spacings a, using a series of configurations sepa- 
rated by 10 supersweeps. In analysing our data we ruled 
out any measured value of the persistence amplitudes 
with relative error larger than 0.5. The persistence am- 
plitude data was fitted to the formula 

{y\e- HT \V) = e- v ^ T \(n\y)\ 2 , (4) 
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FIG. 1: Overlap between the axial state and the ground state, 
at various lattice spacings, plotted against the number of lat- 
tice points n = r/a between charges (for simulation details, 
see the text). The overlap is independent of /3. 









1 1 


1 1 1 1 1 


1 1 1 




- 






P = 2.325 


- 








□ P = 2.375 










O P = 2.425 










A- A p = 2.475 


— 








V P = 2.525 








m * 


[>--|> p — 2.3/D 








<S P = 2.625 




— 






>^< P = 2.675 


— 


- 

_ 




? A 3£ 




- 








i 










\w a. $ 


TIT 








v v 7 "a... 


. h 










A... 
















1 , 


1,1,1 


1,1 





I I I I I I I I I I I I I I I I I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



r[fm] 

FIG. 2: Overlap between the axial state and the ground state, 
plotted against the physical separation of the charges. Some 
lines have been added to guide the eye. It is seen that, for all 
r, the overlap is smaller for finer lattice spacings (larger 0). 



which assumes that the elapsed time T is large enough 
for contributions from excited states to decouple, follow- 
ing ((T|) . The axial state data we present was derived from 
measurements taken at an elapsed time of T > 5a, for 
which no significant deviations from ^ were observed. 
The overlaps were extracted from the logarithm of the 
data by a weighted least-squares fit to a straight line 
(measured for a fixed value of spatial separation r of the 
static charges). The quality of such a fit can be con- 
trolled by inspection of the values of two parameters: Q, 
the goodness-of-fit, and the reduced x 2 value, \ 2 ' jv where 
v = N — 2 denotes the number of degrees of freedom for a 
linear fit of N data points (this is not to be confused with 
the state |x))- O ur acceptance criteria were x 2 / v < 3 
and Q > 1CT 3 . 



B. Numerical results 

In Figure[T]we plot the overlap of the axial and ground 
states against the number of lattice sites n between the 
charges. Only data points for the overlaps with relative 
errors smaller than 0.35 are displayed. We note first that 
all data points lie on top of one another, signalling that 
the overlap is independent of the lattice spacing (and 
therefore 0) used for the simulations. Given that the 
U(l) overlap, ([3]), is a function of rA, we may expect 
that the SU(2) overlap is a function of r/a, as 1/a is the 
UV cutoff on the lattice. We therefore fit the data to an 
exponentially decaying function of the form 

|(n| X >| 2 = Ce X p( -A~V (5) 



The fit is shown as the solid line in Figure Q] We observe 
a very close fit to the data, for the parameter values C — 
1.307, A = 0.359. 

For any given charge separation r, the continuum limit 
a — y corresponds to n — > oo. Our simulations clearly 
show us that as n increases the overlap of the axial state 
with the ground state drops exponentially. Both this 
exponential decay and the dependence on r/a parallel 
the U(l) results. As we approach the continuum limit, 
a — y (which corresponds to /3 — > oo), we probe more of 
the ultraviolet artifacts of the infinitesimally thin string 
which lowers the overlap with the ground state. 

The overlap may also be expressed as a function of 
the physical separation of the charges. The conversion 
was performed by interpolating between known values 
for the string tension aa 2 ((3) at various values of (3. The 
string tensions and the perturbative 1-loop interpolation 
formula were both taken from [2l|. Figure [2] displays the 
overlap as a function of r at various (3. As (3 increases 
the overlap decreases for all separations of the charges. 
If we fix r at a given separation and take the continuum 
limit, (3 — y oo, the overlap again tends to zero. If the 
trends observed in our results continue to hold at smaller 
lattice spacings, then in the continuum limit the overlap 
with the ground state will vanish. 

Recent papers 0, [23| have made two claims regard- 
ing the geometry of the SU(2) flux tube: firstly, that at 
a fixed lattice spacing the width of the flux tube grows 
with increasing charge separation but, secondly, that at 
finer lattice spacings thickening of the tube is a sublead- 
ing effect, as the width of the tube is proportional to the 
lattice spacing. The first claim agrees with the interpre- 
tation of our above results - we observe that the very 
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FIG. 3: Smeared state overlap with the ground state at various lattice spacings. For a fixed smearing level S the results lie on 
the same curve, they are /3 independent. We have added data for the axial state (S = 0) for comparison, and some lines to 
guide the eye. Lattice: 500 configurations, T = 4a, relative error of Wilson line correlator < 0.5, x 2 /u < 3 and Q > 10~ 2 . 



thin string is a poorer description of the ground state at 
larger separations. The second claim would imply that 
in the continuum limit all thickening of the tube is sup- 
pressed, and that the flux tube of SU(2) would in fact be 
an infinitely thin string of flux. However, we have seen 
that for a given r the overlap of this axial state with the 
ground state drops exponentially as we move toward the 
continuum limit. 

The linear a dependence of the flux tube width de- 
scribed in (22l . [23l | was observed to set in at around 
a ~ 0.06 fm, or (3 ~ 2.600. Our largest (3 values 
(/? = 2.675) probe this region, yet we see no change to the 
functional form of exponential decay. We predict that in 
the continuum limit the overlap between the axial and 
ground states vanishes, and that the SU(2) ground state 
is not well described by a thin string stretched between 
the charges. 



III. SMEARED STATES 
A. Smearing the string 

If a thin string is not a good description of the q-q 
ground state of the glue, then we must look for a state 
which describes a more dispersed distribution of glue. 
We have seen in U(l) theory that it is the infinitesimal 
thickness of the string which leads to a vanishing over- 
lap with the ground state, and a similar interpretation is 
supported by our numerical results in SU(2). As thick- 
ening the string will soften the UV behaviour we would 
therefore hope to see an improvement in the overlap with 
the ground state. 

A state with a thicker string may be prepared on the 



lattice using "smearing" [2H HH, [26[ . Smearing replaces 
links by a sum of their adjacent staples projected onto 
an SU(2) group element. The smeared axial state is pre- 
pared by smearing the links between the charges, the ef- 
fect being to broaden the string of glue. It is well known 
that the use of such smeared operators greatly improves 
the accuracy of calculations of, for example, t he g lueball 
spectrum [2 71 ] and the inter-quark potential [28j. This 
is because smearing an operator reduces its sensitivity 
to higher excitations of the theory, improving the pro- 
jection onto the ground state in the calculation of large 
time Wilson loops. This is the origin of the alternative 
name 'overlap enhancement'. We expect, then, that re- 
moving the UV modes from our state by smearing the 
string should give us a state which is closer to the true 
ground state. 



B. Numerical results 

Repeated smearing gives an increasingly smoother and 
more dispersed configuration. We smeared our axial state 
with a number S of smearing steps from S — 1 to S = 10. 
Wilson loops with T > 4a were sufficient to project onto 
the ground state and 500 lattice configurations were used 
in the simulations. 

Our results are plotted in Figure [3l as a function of 
r/a, and in Figure 2] as a function of the separation r. 
These plots also show some axial state data (no smear- 
ing steps, i.e. S = 0) for comparison. For a given sep- 
aration, smearing clearly improves the overlap with the 
ground state. Figure [3] shows us that the overlap is (3 in- 
dependent at fixed S (as it was for the axial state, S = 0). 
For a fixed r, we approach the continuum limit by fol- 
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FIG. 4: Smeared state overlap with the ground state, plot- 
ted against physical separation. Axial state data (S = 0) is 
displayed for comparison and some lines have been added to 
guide the eye. Legend as in Figure [3] 



lowing one of the curves in Figure [3] to the right. Each 
curve falls as we go to the continuum so the overlap drops 
to zero. Although not as severe as the exponential drop 
of the axial (S = 0) state, the overlap at a given r still 
decreases as we go to finer lattices, suggesting that the 
smeared overlap also vanishes in the continuum limit. 

The overlap enhancement provided by a given S is 
therefore fixed and does not increase as we go to finer 
lattices. The practice of smearing to, say, 10 steps, is 
an artifact of our current lattice sizes. As computations 
on finer lattices become more common, a given level of 
smearing will not produce the same improvement in the 
overlap as for current lattices. Instead, the smearing level 
will have to be increased to see a comparable improve- 
ment. 

So far, we have considered states which describe a thin 
string, or a thicker string produced from smearing. We 
have seen that the thin string gives an increasingly worse 
description of the ground state as we go to the continuum 
limit. We have also seen that the amount of overlap 
enhancement provided by smearing depends not on (3 but 
just on the number of smearing steps. The overlap of the 
smeared state also tends to zero in the continuum limit. 

Higher numbers of smearing steps describe successively 
more widely spread gluonic fields, and increasing the 
number of smearing steps increases the overlap with the 
ground state (at a given f3). We will therefore now turn 
to a model which we expect to describe more widely 
distributed glue. We will find a markedly different be- 
haviour of the overlap in the continuum limit. 



IV. THE COULOMB STATE 

A. Construction and properties 

The axial state © describes a heavy meson without 
gauge invariant constituents. Dressings can, though, also 
give us individually gauge invariant colour charges, de- 
scribing the expected short-distance physics of QCD. Wc 
will now examine the dressing which generalises Dirac's 
static electron [3| to a static, heavy quark in non-abelian 
gauge theories |2|. A gauge invariant charge can be cre- 
ated from the vacuum by the operation 

h- 1 [A,x}q(x)\0) , (6) 

where h~ x is a functional of A obeying [3. l29j 



1 

— / 

fJ 



d j [h- s -A j h+-h- l d j h ) = 



(7) 



This equation is to be understood in a Schrodinger repre- 
sentation where we have diagonalised the operator A. It 
defines h[A, x], the field dependent transformation which 
takes A into Coulomb gauge. The dressing itself trans- 
forms under gauge transformations as h~ x — > h~ l U if 
q — > U~ l q, giving a gauge invariant fermion. Equation 
(O may be used to perturbatively construct this Coulom- 
bic dressing by expanding in powers of the coupling, giv- 
ing a gauge invariant charge with a well defined colour 

0, |29] . A dramatic simplification occurs if we choose 
to work in Coulomb gauge, however, as the dressing be- 
comes a trivial factor of unity in colour space @, [29| . 

1. e. 



h-^AtX] 



1 



(8) 



diA i= 



This is easy to see from the definition of h[A, x] as the 
rotation into Coulomb gauge - if we arc already in that 
gauge, no rotation is needed. Our non-abelian, colour 
singlet Coulomb state | $ } contains two such gauge in- 
variant charges separated by a distance r. In Coulomb 
gauge we therefore have 



\*) = q{x 2 )q(x 1 )\0) 



(9) 



traced over colour indices. In perturbation theory, the 
gluonic fields around these charges are distributed over all 
space, and the potential (<I> \H\ <I>) between the charges 
is, at lowest order, 

a 2 C 1 

( $ \H | $ ) = + self energies + 0{g A ) . (10) 

4tt r 

This is the lowest order contribution to the familiar 
interquark potential. At higher orders, screening and 
anti-screening structures emerge |3(ij . It is useful here 
to outline the properties of the corresponding state in 
U(l). There, if we have only heavy charges, the Coulomb 
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state is the ground state in the charge-ant icharge sec- 
tor. It describes two individually gauge invariant charges 
surrounded by Coulomb fields. When we allow light 
fcrmions, however, we know that the effect of virtual 
pairs is to screen the charge, which lowers the energy. 
This results in an additional, but gauge invariant, con- 
tribution to the U(l) Coulomb state which incorporates 
the screening effects [3l| and has the form of Polyakov 
lines. Returning now to the non-abelian theory, it is 
well known [32j that glue both screens and anti-screens 
the heavy charges. The dressing /i _1 [A,a?] contributes 
the anti-screening effects [331 ] - The additional, gauge in- 
variant contributions which give screening effects come 
from short Polyakov lines generated as the state evolves 
in time. 

Our Coulomb dressing describes individually gauge in- 
variant charges in perturbation theory, where we can de- 
scribe the physics of our dressings analytically. Non- 
perturbatively we know that any description of a physical 
colour charge must break down, as no such objects are 
observed (as asymptotic states) in nature. The Gribov 
ambiguity has been shown to generate just this break- 
down 0, [29| . Confinement will appear at the scale where 
Gribov copies reintroduce an unphysical gauge depen- 
dence to the single charge dressing. Beyond this scale 
only dressings for hadronic (colourless) states exist. This 
could force our two individual dressings, h[A, a^] for 
q(x2) and /i _1 [A,a;i] for q(x\), to combine into a single 
mesonic dressing, i.e., the Coulomb state could still, non- 
perturbatively, confine the charges 0, In fact, it is 
known from lattice results that the non-abelian Coulomb 
gauge potential, which is an upper bound to the full po- 
tential, is confining with a linear rise [l3l . |34| , in contrast 
to the abelian theory. 

We now turn to the numerical preparation of our state 
and the calculation of its overlap with the ground state. 
We will see that, non-perturbatively, the Coulomb state 
is a good approximation to the true ground state, and 
that their overlap in fact improves as the lattice regulator 
is removed. 



B. The Coulomb state on the lattice 

Our simulations were performed using an improved 
lattice action |2l[ which gives good scaling and rota- 
tional symmetry properties. Calculating the persistence 
amplitude for the Coulomb state requires bringing the 
lattice into Coulomb gauge and evaluating correlators 
of Polyakov lines in the gauge fixed configuration [lj|. 
We used an iteration-overrelaxation algorithm to fix the 
Coulomb gauge by maximising the gauge fixing func- 
tional 
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FIG. 5: Improved action simulation of the Coulomb overlap 
at fixed (3 = 1.55, 20 4 lattice. As T increases contributions 
from excited states decouple, and we see a decrease in the 
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plotted - however, the difference between results at different 
T is decreasing. 



During gauge fixing we monitored the behaviour of 

N 3 

4.Y 



1 N 3 



(12) 



x=l b=l 



Here B(x) = Y?i=A u i( x ) + U H X - and r b arc the 
three generators of SU(2). Expanding in the lattice spac- 
ing we find 



ti[T b B(x)} = a 2 diA\{x) + 0(a 3 ) , 



(13) 



which shows that S 2 is a measure of violation of the 
Coulomb gauge condition. As a stopping criterion we 
demanded 6 2 < 10~ 10 . Our simulations cover quark sep- 
arations from roughly 1.8 fm down to 0.05 fm. 

In a first run, we studied the T behaviour of the data 
using a 20 4 lattice and (3 = 1.55, with the aim of identify- 
ing T min such that excited states effectively decouple for 
T > Train- The results of this run are shown in Figure O 
which plots the overlap as a function of r for several val- 
ues of T min . 

Although the static potential is stable for T m in > 3, the 
overlap retains a T m i n dependence. The high value of the 
overlap, over 70%, and the sensitivity to T mm indicate 
that the Coulomb state has a significant overlap not only 
with the ground state but also with low lying excited 
states, which have not decoupled. We must therefore 
go to larger T values to see the overlap with only the 
ground state. It seems, though, that the higher modes 
will decouple not far beyond T = 5a, as our results show 
that the difference between overlaps decreases as T is 
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FIG. 6: Improved action simulation of the overlap between the 
Coulomb and ground states. At a fixed time, here T = 3a, 
we observe an increase in the overlap in the continuum limit 
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FIG. 8: Wilson action overlap of the Coulomb state and 
ground state, plotted against the number of lattice sites 
n — r/a between the charges. For lattice data, see the text. 
Measurements were taken from T > 3a. Notably, the overlap 
increases as we go to the continuum, (3 — > oo. 
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channel. 

We plot in Figure [S] the overlap at various lattice spac- 
ings and fixed T. Figure [7] summarises the data, showing 
both the T and j3 dependencies. The difference between 
the behaviour of the Coulomb and axial states in the 
continuum limit is striking - for a given r the Coulomb 
overlap increases with decreasing lattice spacing, becom- 
ing a better description of the true ground state as we 
approach the continuum. Does this pattern hold as we 
go to finer lattices? To answer this question and study 
our Coulomb state closer to the continuum limit we now 
turn to the Wilson action, which will allow us to study 
the overlap at larger f3 values. This will also show that 
our results are independent of the discretisation method 
employed. We present the results of these simulations 
below. 



Wilson action results 



FIG. 7: A summary of the improved action results for the 
Coulomb overlap, plotted at a range of (3 and T values. 



increased, and the curves at T = 4a and T — 5a are 
already very close together. The results at our larger T 
values should therefore serve as a good approximation 
to the true overlap. Having addressed the residual In- 
dependence, we now study the overlap for varying (3. We 
note first, however, that the dependence on T m - ln suggests 
a small or possibly zero mass gap in the quark-antiquark 



Our Wilson action results for the Coulomb overlap 
are shown in Figure [5] as a function of r/a. We ob- 
serve a markedly different behaviour to both the axial 
and smeared states. Firstly, we note that the results do 
not all lie on the same curve. The form of the data no 
longer fits a (3 independent function of r/a, as found for 
the axial overlap, Figure [TJ Instead, as (3 increases the 
Coulomb overlap increases for all values of r/a. The ori- 
gin of the axial overlap's dependence only on r/a was the 
UV physics of the infinitesimally thin string of glue be- 
tween the charges. The configuration of glue described 
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FIG. 9: Wilson action overlap of the Coulomb state and 
ground state, plotted against the physical separation of the 
charges. The overlap increases as we increase (5. 



by the Coulomb dressing has no such structure. This 
eliminates the UV artifacts of the axial state and so the 
Coulomb overlap depends on both r and a separately 
rather than on only their ratio. 

We now turn to the overlap as a function of the phys- 
ical separation of the charges, plotted in Figure El anal- 
ogously to Figure [S] Notice that the order of the curves 
corresponding to simulations at different lattice spacings 
is reversed as compared to the axial plots. At short sep- 
arations the overlap with the ground state is almost per- 
fect and drops more slowly with increasing r than the 
exponential drop of the string state (see Figure [2]) . It 
is also almost independent of the lattice spacing in this 
regime. As we draw the two charges apart, our physi- 
cal expectation is that a flux tube forms, confining the 
charges into a meson. In this situation the Coulombic 
description may not be expected to be appropriate, and 
indeed we see that the overlap drops with increasing dis- 
tance. However, this trend does not continue indefinitely. 
For a given lattice spacing we see, at some physical sep- 
aration, a levelling off in the overlap. This phenomenon 
occurs at different distances for different lattice spacings 
- as the lattice spacing is decreased, and we move to the 
continuum limit, this levelling of the overlap begins at 
smaller and smaller physical separation, giving a good 
and almost r independent overlap with | fi ) . The Wil- 
son action results therefore confirm our previous findings 
with the improved action, Figure [71 

The most significant difference between the axial and 
Coulombic states is that the overlap of the former (and 
its smeared counterparts) decreases as we go to the con- 
tinuum limit, for all separations of the charges. The 



Coulombic state, however, has a good overlap with the 
ground state which increases in the continuum limit. For 
any charge separation r, we move up the graph in Figure 
[5] away from and towards 1. 

The small difference from 1 in the overlap of the 
Coulomb and ground states is reminiscent of perturba- 
tive results in (Coulomb gauge) time-independent per- 
turbation theory [32], l35| . The overlap between the g = 
ground state (essentially AT? — 1 copies of the abelian the- 
ory) and the perturbed q-q sector ground state goes like 
1 — 0(g 2 ) for g small. Here we have seen that the over- 
lap between our Coulomb state and the true ground state 
differs, non-perturbatively, by a small quantity which de- 
creases as we go to the continuum. We therefore expect 
that the physics of the Coulomb state should give a good 
description of the true ground state physics, implying in 
particular that the ground state of the glue around heavy 
charges takes a form much thicker than that of a thin 
string. Due to non-perturbative effects, such as the Gri- 
bov ambiguity, we cannot rule out the possibility that, 
non-perturbatively, the Coulomb state becomes the true 
ground state in the continuum limit - as discussed, the 
non-abelian Coulomb potential is linearly rising for large 
separations and a good agreement between the Coulomb 
string tension and the full string tension was reported 



V. CONCLUSIONS 

We have constructed and analysed models of the q- 
q sector ground state in Yang-Mills theory. Our states 
are gauge invariant and describe different distributions 
of glue around the fermions. 

Our first ansatz, the axial state, modelled the flux tube 
by a string of glue stretched between the fermions. We 
calculated the overlap between this state and the ground 
state in SU(2). We saw that for any separation of the 
fermions the overlap decayed exponentially with the ul- 
traviolet cutoff provided by the inverse lattice spacing. 
These results parallel those of U(I), where the infinitesi- 
mal transverse extension of the string leads to a vanishing 
overlap between the axial state and the ground state in 
the continuum. We conclude that the SU(2) flux tube is 
not well described by a very thin string. 

This should be contrasted to the confining phase in 
compact U(l) where, at a fixed lattice spacing, we have 
seen [2(| that the axial state provides a good description 
of the ground state. Investigating whether this statement 
holds in the continuum limit could shed light on the dif- 
ferences between abelian and non-abelian confinement 
(see also (36J). It is also interesting to speculate on the 
role of the SU(AT C ) axial state for large N c . Here Yang- 
Mills theory is expected [13] to have a dual description 
as a theory of strings, and it may be that in such a limit 
the axial state again becomes a better description of the 
ground state. 

Smeared states are formed on the lattice by replacing 
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links with sums over staples. The effect of this procedure 
is to broaden the string, removing the UV modes. We 
found that the overlap enhancement offered by smear- 
ing is independent of [3 for a given number of smearing 
steps. Our results show that, for a given separation of 
the charges, the overlap between the smeared state and 
the ground state still drops toward zero in the contin- 
uum limit, though it does so more slowly than for the 
unsmeared axial state. This implies that for larger lat- 
tices more and more smearing will be required to improve 
simulations. 

Our final model was the Coulomb state, which de- 
scribes, in perturbation theory, two individually gauge 
invariant colour charges surrounded by non-abelian 
Coulombic fields. This ansatz has through its relation- 
ship to Coulomb gauge fixing 0, Hsj and its Gribov 
copies, a rich non-perturbative structure. We observed 
that it has a much better overlap with the ground state 
than the axial descriptions did, which implies it captures 
more of the true ground state physics. Since the overlap 
actually increases as the lattice spacing is reduced, we 
cannot exclude the possibility that the Coulomb state 
is indeed the true ground state of the q—q sector in the 



continuum limit. 

There exist many further questions to explore - for ex- 
ample, it would be interesting to study the role of Gribov 
copies in our numerical calculations and their contribu- 
tions as a function of the lattice spacing. The result that 
the Coulombic description of the ground state yielded 
the best overlap at small lattice spacings deserves future 
study, using large scale simulations to further approach 
the continuum. Our understanding of the geometry of 
the flux tube would be increased by constructing trial 
states with a transverse profile for the tube and max- 
imising their overlap with the ground state. It would be 
also be very interesting to explore our surprising result 
that the data supports a small, or zero, mass gap in the 
meson sector. 
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